Can satellites deliver substructures and black holes to inner halo 

by dynamical friction? 
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ABSTRACT 

A fully analytical formulation is developed to make dynamical friction mod- 
eling more realistic. The rate for a satellite to decay its orbit in a host galaxy 
halo is often severely overestimated when applying the ChandraSekhar's formula 



lO , without correcting for the tidal loss of the satellite and the adiabactic growth of 

the host galaxy potential over the Hubble time. As a satellite decays to the inner 



and denser region of the host galaxy, the high ambient density speeds up their 
i exchange of energy and angular momentum, but shrinks the Roche lobe of the 

satellite by tides. Eventually both processes of orbital decay and tidal stripping 
hang up altogether once the satellite is light enough. These competing processes 
can be modeled analytically for a satellite if we parametrize the massloss history 
by an empirical formula. We also incorporate the adiabatic contraction of orbits 
j§ ! due to growth of the host potential well. Observed dwarf galaxies often show a 

finite density core, which determines how much inwards its remnants could be 
^ \ delivered to the host galaxy. Ghost streams or remnant cores or globular clusters 

^ ! of satellites should populate preferentially the outer halo (e.g., the Magellanic 

stream and the Ursa Minor dwarf), rather than the inner halo (e.g., the Sgr 
stream and the Omega-cen cluster). Massloss due to strong tides in the inner 
galaxy also makes it problematic for any small central black holes in nuclei of 
satellite galaxies to decay in orbit and merge into the center of the host galaxy. 



1. Introduction 

Current theory of galaxy formation favors the idea that galaxies form hierarchically 
by merging smaller lumps or satellites. Dynamical friction dissipates the orbital energy 
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of the satellites so that they sink deep into the host galaxy potential well, where they are 
disintegrated and virialized via baryonic feedbacks and tidal stripping. These processes might 
have determined the density profile of virialized halo of the host galaxy (Syer & White 1998, 
Dekel et al. 2003). 

There are about 150 and 300 globular clusters, and a few dozen dwarf sateUites of mass 
1O^~^M0 in the Milky Way and M31 respectively. It is tempting to associate dwarfs satellites 
and globular clusters in the Milky Way as the markers/remnants of past hierarchical merging 
events. Indeed there are several examples of possible streams of remnants in the Milky Way 
(Lyndcn-Bcll & Lyndcn-Bell 1995) including the recently found Galactic ring or Carnis 
Major dwarf galaxy, traced by a grouping of globular clusters (Martin ct al. 2003). A giant 
stream is found in the Andromeda galaxy (McConnichie et al. 2003, Ferguson et al, 2002). 
Among these the Sagittarius dwarf galaxy stream (Ibata et al. 1997) between radius 10-50 
kpc from the Galactic center is perhaps the best example. It brings in at least 5 globular 
clusters to the inner halo, including M54, the 2nd most massive cluster of the Milky Way. A 
still mysterious object is a;-Gentauri , with about lO^L© at 5 kpc from the Galactic center 
which has the morphology of a globular cluster, but has multiple epoches of star formation 
and chemical enrichment (see Gnedin et al. 2002 and the c<j-Centauri symposium). Another 
example is the Gl cluster, a very massive globular-like object with about IO^Lq at about 40 
kpc from the M31 center (Meylan et al. 2001). A system (NGC1023-13) almost identical to 
Gl is also found in the SO galaxy NGC1023, at a projected distance of about 40 kpc from the 
host galaxy (Larsen 2001). Freeman (1993) suggested that such systems are the remnants of 
nucleated dwarf satellites, with their outer tenuous dark matter and stars being removed by 
galaxy tides. M31 also has an unusual collection of clusters as luminous as a;-Centauri within 
5 kpc in projection from its double-peaked center. Given the above evidences or signs for 
infalling objects in our galaxy and M31, it is interesting to ask whether some of the inner 
globulars could have also been the result of mergers. Beyond the Local Group, minor mergers 
are sometimes speculated as the mechanism to deliver massive black holes and gas material 
into the nucleus of an AGN to account for the directions of the jets and nuclear dusty disks. 
A very interesting related issue is whether giant black holes acquire part of the mass by 
merging the smaller black holes in the nuclei of infalling sateUites. 

For this paper we revisit the basic theoretical questions: what is the condition for a 
dwarf galaxy to decay into the inner halo? What are the possible outcomes of tidal stripping 
of a dwarf sateUite? How often do we get a system like a;-Centauri or a naked black hole 
near the host center? The answers to these questions will help us to test the validity of 
the theory of hierarchical merger formation of galaxies. The key mechanism for satellites 
to enter the inner galaxy is dynamical friction, where the gravity of the satellite creates a 
wake of overdensities in the particle distribution of the host galaxy, which in turn drags 
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the motion of the satellite with a force proportional to m{tY, where m{t) is the mass of 
the satellite. Another process is tidal disruption, where the object sheds mass with each 
pericenter passage, and the remnants are littered along the orbit of the satellite. The above 
two processes compete with and regulate each other: orbital decay increases the tidal field, 
which reduces the mass of the satellite, hence slows down the orbital decay. Some examples 
of these effects have been shown in Zhao (2002) in the case of a;-Centauri . 

The analytical formula of ChandraSekhar is widely used for gaining insights on dy- 
namical friction because of the time-consuming nature the more rigorous N-body numerical 
simulation approach. It is a customary practice in previous works to model the orbital decay 
of a satellite as a point mass of a fixed mass. However, the fixed mass approximation is un- 
necessary and could seriously overestimate dynamical friction because of neglecting massloss 
It is essential in calculations of satellite orbits to model the dynamical friction and 
massloss together since they regulate each other. 

In the past the massloss and the orbital decay are often modeled in the ab. initial 
fashion, resulting coupled non-linear equations without simple analytical solutions. In such 
models, the sateUite mass is often modeled as a function of the sateUite's tidal radius, hence 
various factors come in, including the orbital position of the satellite, the density profile of 
the satellite (e.g., Jiang & Binney 2000, Zhao 2002, Mouri & Taniguchi 2003, Kendall et 
al. 2003). However, these complications are not always necessary since the massloss history 
is rather similar in simulations with very different initial conditions (the mass is generally 
a stair-case like a function of the time), so could be parametrized in an empirical fashion, 
by-passing the uncertain assumptions of the satellite initial profile. This could be useful for 
exploring a large parameter space of the satellite initial condition. 

Another unnecessary approximation but common practice is to use a static potential 
for the host galaxy. This again is unphysical since galaxy halos do grow in hierarchical 
formation scenario partly because of galaxy merging, and partly because of the adiabatic 
contraction of the baryonic disk and bulge; galaxy rotation curves Vdrir-^t) can change by 
significantly before and after the formation of baryonic disks and bulges in mass models for 
the Milky Way and M31 (Klypin, Zhao & Somerville 2002) and in generic CDM simulations 
with baryons (e.g. Wright 2003). The growing gravitational force tends to restrain the radial 
excursions of the satellite while preserving the angular momentum. The dynamical friction 
or drag force is also proportional to the growing density p(r, t) of ambient stars and dark 
particles in the host galaxy. 

In fact, it is conceptually simple to incorporate massloss and growing potential while 
keeping the problem analytically tractable: the deceleration — ^ is simply proportional to the 
satellite mass ^{t) and the ambient density p(r, t) of the host galaxy at the time t. Without 
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massloss ChandraSekhar's formula would predict very efficient braking of the orbits, enough 
to make a high-mass satellite of ~ IO^^Mq (the mass of the LMC or M33 sized object) to 
decay from a circular orbit at ~ lOOkpc to the very center of a high brightness galaxy in a 
Hubble time, delivering remnants into the inner galaxy. Here we study the effect of massloss 
and growing potential on the result of orbital decay, and the distribution of remnants. We 
present fully analytical results for calculating the decay rate for satellites on eccentric orbits 
in a scale-free growing isothermal potential. 

The structure of the paper is following: in S2 we give the analytical formulation of the 
problem, in S3 and S4 we present results of application of our analytical model to globulars 
and dwarf satelhtes, in S5 we discuss and in S6 we summarize. 



2. Analytical model of orbital decay of a shrinking satellite in a growing host 

halo 

Consider a satellite moving with a velocity Vg on a rosette-like orbit in a spherical host 
galaxy potential 0(r, t) with a rotation curve Vdi- Assume it has an initial mass rrii, the 
orbital decay from an initial time t = ti to the present day t = to can be modeled using 
ChandraSekhar's dynamical friction formula (Binney & Tremaine 1987) 



where tftc is the instantaneous dynamical friction time. Manipulating the equation, we find 
the (specific) angular momentum j(t) of the satellite decays as 

Here the dynamical friction time is given by 

t- = [47rGp(r,t)]^e, (3) 

cir 

where m{t) is the mass of the sateUite, p(r, t) is the density of host galaxy at the orbital 
radius r{t), and ^ is some dimensionless function of the rescaled satellite speed ^y^^- 

Here we adopt a Singular Isothermal Spherical (SIS) host galaxy model with a rotation 
curve generally growing with time Kir(i); it is normalized by the present time value Kir(io) = 
Vq. In this model the potential and density are given by 

Hr,t) = Vl{t)lnr, p(r,t) = ^. (4) 
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Since the stars and dark matter velocity distribution in the host galaxy is an isotropic 
Gaussian, the dimensionless velocity function 

2u 

erf (m) — — ^ exp {—u^) 



In A, u{t) ^ ( ^ ) , (5) 



an) 

where the dimensionless Coulomb logarithm In A is typically between unity and ten. 



Combining the above equations we find the decay rate for the angular momentum is 
given by 

dj(t) j{t) 1 G'/i(t)Kir(i) 



dt tfrc ' ifrc j{t)'^ 

where we define an effective bound mass of the satellite at time t 



(6) 



li{t) = m{t)i{uy cos2(a), u{t) = ^y^, (7) 

where u{t) and a{t) are the rescaled speed and the pitch angle of the orbit at time t. As we 
can see, the effective mass lumps together several time-varying factors, including the 
satellite massloss and orbital epicyles. Introducing this effective mass simplifies the physics 
details because the braking rate of specific angular momentum —dj{t)/dt is now simply 
proportional to 

Integrating over time, we find the change in the specific angular momentum between 
time ti to to is given by 

-Af = fit,) - fito) = 2Gm,VoT, r ^ {t, - t^h ^ ^ ^^^^^dt, (8) 

where we define 7 as the reduction factor from fixed satellite model, and define r as the 
effective duration of dynamical friction, m,i and Vq are the mass at initial time t, and the 
rotation curve at the present to- The above formula allows us to calculate the evolution of 
the specific angular momentum j{t) of a generally eccentric satellite orbit. 

It is interesting to ask what the necessary condition is for an outer halo satellite to reach 
the inner galaxy or even the center. These satellites will be destroyed with their remnants 
in the inner galaxy, e.g., by the disk and/or bulge shocking. To facilitate the comparison 
with observations, it is sometimes better to define an instantaneous orbital size ^ 

S{t) ^ ^-P- (9) 



^We can make this defination irrespective of the eccentricity of the orbit, which we does not enter our 
calculation explicitly. Roughly speaking the orbital size S{t) is the geometrical mean of the pericenter radius 
and apocenter radius in a static potential. 
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at time t. Let the present orbital size So — ^ < i?disk ~ 15kpc, and define j^isk — -Rdisk x Vo 
as the critical angular momentum to reach within the disk at the present time. Manipulating 
eq. (8) we find the initial angular momentum ji or orbital distance Si — ^ must satisfy 

0<So< i?disk, r ^ r ^^^^^^f^dt. (10) 

In other words, satellites initially on orbits larger than Si would not be able to deliver a 
globular or dwarf galaxy to the inner galaxy irrespective of mass loss assumed. 

Eq. (8) is our main analytical result. There is also an interesting geometrical interpre- 
tation to it. Consider a satellite of mass m{t) on a circular orbit in a static SIS potential, 
we have 

lnAf«2.5, n{t) ^m{t)^{l) ^QA2hiAm{t) ^m{t), (11) 

where we took a typical value for the Coulomb logarithm (e.g., Panarrubia, Kroupa & Boily 
2002). We can rewrite the angular momentum equation (eq. 8) as 

o, 2nG(m) ^l'' m{t)dt 

IttA^^I « (FoAt) , (m) = \ (12) 

to — li 

where the l.h.s. is the area swepted by the decaying orbit, and in the r.h.s. 2tiG (m) V^'^ 
is the "circumference of influence" of the satellite, and the Vq/^I = Vo(to ~ ti) is the length 
of the satellite's orbital path, and the multiplication of the two is the area swept by the 
circumference of influence of the satellite. The above equation implies that the two areas are 
comparable. Interestingly the decay of the orbits depends on the sateUite mass through the 
approximation AS^ oc (m) At, i.e., it is the average mass of the satellite that determines the 
rate of orbital decay. Finally note eq.(8) applies to the evolution of specific angular momen- 
tum (S.A.M.) of an eccentric orbit in a time- varying potential. A complete description of the 
orbits should also include the evolution of the orbital energy, which unfortunately is more 
complex analytically, and is not studied in detail here. It is not essential for our conclusion, 
but is perhaps convicnent to assume efficient dynamical decay of the orbital energy during 
the pericentric passages hence the radial motion is damped and orbit circularizes at the end. 
Adiabatic growth of the potential also tends to circularize the orbits. 



Si^ 



2GmiT 



2.1. Models for a growing host potential and a diminishing satellite 



Classical singular isothermal model assumes a fixed potential or a time-invariant rotation 
curve. In reality galaxies grow substantially from redshift of a few to now. The growth of 



- 7- 



the dark halo is scale-free in the hierarchical scenario, so it is plausible to approximate the 
growth of the rotation curve as a power-law in time as follows, 



Kir(^) 



0<t<tn 



(13) 



where the rotation speed at some earlier time < t < to is generally smaller than the 
present day (t = to = 14Gyr) value Vq. At the present epoch Vcu{to) = Vq = 200kms~^ is 
appropriate for the Milky Way. A few sample evolution models are shown in Fig.l. A value 
of I < p < I imphes an evolution of 10%-25% in rotation speed from redshift one to now, 
which seems reasonable in models of the Milky Way before and after the formation of the 
disk and the bulge (Klypin, Zhao, Somerville 2001; Wright 2003). 

Previously the tidal stripping of the sateUite is either modeled by N-body simulations 
(e.g., Tsuchiya et al. 2003, Johnston et al. 1999) or semi-analytically using tidal radius 
criteria (e.g., Zhao 2002, Kendall, Magorrian & Pringle 2003). ^ In both approaches one 
needs to assume a rigorous description of the mass profile of the progenitor satellite. It is 
unclear whether this approach can model efficiently the realistic large scatter in the circular 
velocity curves of observed dwarf galaxies (see, e.g., Fig.4a). 

Here we take a very different approach. We model the mass of the progenitor galaxy as 
a simple function of time with two free parameters. We do not need an explicit prescription 
of the satellite density profile. Motivated by the massloss history typically seen in N-body 
simulations, the satellite mass is modeled to decay in a rate between exponential and linear 
massloss. To simplify the calculations, we lump together several uncertain factors, and 
simply assume the effective mass of the satellite is shed following the following recipe. 



sin2iV7rt\ 



rrin 



2N7r 



ml 



t-t, 
to — ti 



(14) 



where mo and mi are the mass at the present t = to 
in t = ti, and the quantities < t < 1 and N ~ ; 



14Gyrs and when the satellite falls 
100 are the rescaled dimensionless 
time and the number of pericentric passages between time to ^-nd tj. The parameter n is 
tunable with the n — 1 model having a constant rate of massloss, and the n — > oo having an 
exponential massloss; and n = — 1 is a roughly power- law decay. In simulations we typically 
see somewhere in between these cases: < n < 1 for Plummer or King model satellites with 



"^Recently the latter approach has also been extended to model the orbital decay and evaporation-induced 
massloss in the dense star cluster near the Galactic center (Mouri & Taniguchi 2003, McMillan & Portegies- 
Zwart 2003). 
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a sharp fall off after an initially linear massloss (e.g., Panarrubia et al. 2002); for isothermal 
satellite models, the massloss is close to linear n = 1 (Zhao 2002). Note for a completely 
disrupted satellite mo ~ 0. In N-body simulations satellite losses mass mainly in bursts 
near pericentric passages, so the mass is a descending staircase like function of time. This 
is simulated fairly well by our formula. A few sample massloss histories are shown in Fig.l 
where we assume iV = 10 pericentric passages from ti = 4Gyr to now to = 14Gyr. 

Dynamical friction can be greatly reduced by satellite massloss and host growth. The 
reduction factor 7 is computed by substituting eqs. (13-14) into eq. (8). For a rough 
estimation of the reduction factor 7, consider a Gedanken experiment where a satellite enters 
a growing host galaxy at time ti = right after the big bang and is completely dissolved 
(mo = 0) by the time to. Assume the orbits are circularized, so N ^ 00. We find 7 = 
(p+l+i)! = I satellite losses mass hnearly with time (n = 1) in a linearly growing halo 

(p = 1). This estimate is perhaps to the extreme. In reality the formation and mergers 
of the satellites are probably over an extended period of the Hubble time, perhaps starting 
around redshift of 1.5 {ti = 4Gyr), ending around now (to = 14Gyr). Hereafter we consider 
mostly models with tj — 4Gyr, and to — 14Gyr, ^ < p < | and mo "C m^. 



3. Results of application to globulars and dwarf satellites 

Section 2 gives the formalism to predict analytically the evolution of the angular mo- 
mentum of the satellite for a satellite with any massloss history on an eccentric orbit around 
a time- varying potential. 

Clearly not all sateUites could reach the galactic center as a globular or a naked massive 
black hole. Dynamical friction is basicly turned off if the sateUite bound mass drops below 

IO^Mq before reaching the inner galaxy. To reach the inner, say, 15 kpc of the host galaxy, 
which is roughly the truncation radius of the outer disk of the Milky Way, a satellite must 
have presently a specific angular momentum 

< io < jdisk = i?disk1^o = 15kpc X 200 km s-\ (15) 

In comparison the present specific angular momentum of some of the known satellites are 
given in Table 1. The Magellanic stream and Ursa Minor have jo ~ 15000 km s~^kpc much 
larger than jdisk- 

Suppose there was a population of hypothetical dwarf satellites in the outer halo with a 
specific angular momentum comparable to that of the orbits of LMC, Ursa Minor and Fornax. 
Suppose their initial effective mass mj is between Ursa Minor and the LMC (1O^-^°M0) 
at sometime tj between 4 — 14Gyr. We integrate forward in time to answer the question 
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where their remnants are. Clearly the effect of orbital decay is maximized if we take the 
longest evolution time (lOGyr), the highest initial satellite mass and smallest initial orbit 
ji = 60kpc X 250kms~^. This is the case shown by the hatched region in Fig. 2 with the 
satellite starting from the upper right corner and ending to the lower left with a mass IO^Mq. 
The vertical axis is a characteristic orbital distance of the satellite, expressed in terms of 
the SAM (specific angular momentum) divided by a characteristic velocity 200kms~^. This 
orbital distance is roughly the geometrical mean of the apocenter and pericenter. Models 
are shown in order of increasing dynamical braking. The upper shaded zone are models 
with between exponential and linear massloss and a moderate evolution of the potential 
(oo > n > l,p — the lower shaded zone are models with between linear to accelerated 
massloss and a minor evolution of the potential (l>n>0.3,p=|). Qualitatively speaking 
the orbital decay appears to be only modest in all cases, the remnants are generally not 
delivered to the inner galaxy. 

The condition for a satellite to deliver a low-mass substructure to the inner halo (cf. 
eq.l5) or a IO^Mq black hole to the galaxy center is summarized in Fig. 3. The satellite must 
be within 20kpc for the past Hubble time for a low-mass (lO^^^M©) progenitor. It could 
be at a modest distance of 40-50kpc if the progenitor was very massive (10^°Mq) with a 
linear or accelerated massloss (n < 1) and little evolution of the galactic potential (p -C 1). 
Progenitors of the inner halo substructures or central black hole cannot be on orbits of 
specific angular momentum comparable to the LMC, Ursa Minor or Fornax. This illustrates 
the difficulty of making systems such as a;-Centauri as the nucleus of a stripped-off dwarf 
galaxy starting from the very outer halo. Likewise it is difficult for a minor merger to bring 
in a million solar mass black hole to the host galaxy center. The progenitor's orbit must be 
radial and well-aimed at the galaxy center such that the progenitor's tangential velocity is 
< lkms~^ if the satellite comes from an initial distance of IMpc. 

4. Tidal constraints on internal density profile and circular velocity curve of 

satellite 

The basic feature of our dynamical friction models is that they bypass any information 
of the internal density profile of the satellite by specifying m(t) directly. Nevertheless we 
can infer the internal mass distribution, or internal circular velocity curve, of the satellite 
using the tidal radius criteria. More specificly the internal circular motion u^r at the tidal 
radius rt should be in resonance with the satellite's orbital velocity V{t) at the pericenter 
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r{t), meaning that their angular frequencies are equal with 

^-t-J^ (16) 

where tcr is the crossing time, the factor k is the boosting factor of the pericenter velocity 
due to eccentricity, and the velocities 



^M). .v;. = f. (17) 

V n r{t) 

We fix /c = v^3) hereafter, appropriate for a very eccentric orbit with peri-to-apo radius 
ratio of 1 : 5, e.g., the orbit of the Sgr dwarf galaxy; we neglect any circulation of the orbit. 
Use substitutions to eliminate r{t), we have 

"cir J.-1 cir /iq\ 



Gm{t) j{t) ' 

and eliminating Udr we have 

(G«)'oc(^)'=<-' (19) 

where pt is the sateUite's overall mean density at tidal radius r^. This means that as the 
satellite sinks in to the halo with a shrinking mass m{t) and specific angular momentum j{t), 
the satellite's mean density pt increases with the ambient density, hence the tidal radius 
shrinks with the satellite's mass m{t), like what happens with peeling off an onion. The 
tidal peehng-off process effectively maps out the internal mass radial profile of the satellite, 
i.e., m{rt) implicitly. 

Prom the mass radial profile m(rt) we can derive an internal mean density profile pt{rt) 
or more usefully an internal circular velocity profile Udrirt), which is shown in Fig. 3. Here 
we consider only satellites sink in a moderately growing (p = 1/3) host potential, starting 
with a mass and angular momentum similar to the present values of the LMC Depending on 
the assumed massloss history m{t), we have a range of implied circular velocity curves Udr vs 
Tf. At the lower end we have models with rapid (n ^ 1) or exponential {n oo) massloss, 
which show an overall solid-body circular velocity due to an overall uniform density profile; 
such satellites are disrupted rapidly once the ambient density reaches the uniform internal 
density. At the higher end we have models with n <^ 1, which show a nearly Keplerian 
velocity at large radii due to a highly stratified density profile, which behaves like a point- 
mass at large radii; such concentrated satellites are able to retain most of their mass as they 
spirial deep into the denser and denser region of the host galaxy before eventual disruption. 
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The circular velocity curves of all models show a sohd-body core due to a finite central 
density; the faster the rise of the solid-body core, the denser the satellite central density, 
and the further inwards the final disruption location of the satellite. In the middle we have 
models with linear massloss {n = 1), which show a nearly flat velocity curve at large radius 
and a solid-body core. It can be understood analytically why singular isothermal satellite 
losses its mass linearly with time in a singular isothermal host halo of fixed Vdr'- From the 
tidal criteria (eq. 18) and the orbital evolution (cq. 6) we have 

hence we have 

^ oc ^ = est, (21) 
dt dt ^ ^ 

meaning a linear massloss. 

The above formulation is based on the simplifying static picture that the satellite's 
mass is peeled off in successive layers at the shrinking tidal radius. The picture in N-body 
simulations is more complicated, since satellite particles at all radii, e.g., the center of the 
satellite, could be escape at any time. So in principle eq. (19) sets only a lower limit on 
the initial mass profile of the satellite. We shall proceed without making this correction, 
although it is possible to introduce an empirical factor to correct this, e.g., as in Jian & 
Binney (2000). 

Consider satellites evolving from an initial mass 10^°Mq at time ti — 4Gyr to a final 
remnant mass lO^M© at time to — 14Gyr. The mass profiles of the satellites can be predicted 
depending on the massloss parameter n. Several circular velocity curves are shown in Fig.4b 
and Fig. 5. Compare these with the rotation curves of several nearby dwarf galaxies of virial 
mass ~ IQ-'^^Mq (Fig.4a). Observed dwarf galaxies typically have a solid-body core instead of 
a cold dark matter cusp. It is remarkable that our predicted circular velocity curves resemble 
very well those of observed dwarf galaxies, for a variety of initial and final parameters of 
the satellites and the host galaxy. The fact that our models give reasonable mass profiles of 
satellites justifies our empirical parametrization of the massloss history (eq. 14). Generally 
models with n — 1 resemble a cored isothermal model. Models with n » 1 give a solid-body 
circular velocity curve, and models with n <^ 1 give a Keplerian-like curve with a dense 
solid-body core. 

If the parameters of the LMC are typical for sateUites, then from the initial specific 
angular momentum ji ~ Jlmc ~ 50kpc x 250kms~^, we find the initial tidal radius rj is 
around 10 kpc (Fig.4b). The remnant is dehvered to between 30-60 kpc depending on n if 
we also fix p = 1/3. 
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In order to deliver a final remnant to the inner 15 kpc, we must start with an orbit with 
ji < JlmC) hence the initial tidal radius of the satellite Vt < lOkpc. Fig.5 shows the range 
of predictions by varying the parameters p and n. The predicted circular velocity curves 
resemble those of the observed dwarf galaxies (Fig.4a), but have typically high amplitudes 
with typically dense cores. Observed dwarf galaxies have cores of different densities (Fig. 4a). 
The denser is the core of a dwarf galaxy, the further inward is the final disrupted stream, but 
the remnants are never delivered to the very center. Our models suggest that the remnants 
of most of dwarf galaxies should be typically outside the inner 15 kpc of the host galaxy 
halo. 



5. Discussions 

5.1. Effects of Disk and bulge 

There are two limitations of the current analysis. First we assume the effective mass 
is some simple function of time (cf. eq. 15), whereas in reality several complex factors go 
into the modeling of m(t), the velocity function ^{u) and the Coulomb logarithm. Second 
we assume an isothermal dark matter plus stars model throughout the galaxy, hence the 
dynamical friction effect of the disk and bulge are not modeled accurately. However, these 
limitations should not change our qualitative conclusion that the remnants of outer halo 
satellites should remain in the outer halo if tidal massloss is severe. For very eccentric orbits 
the factors in ^{t) may be modulated periodically with the radius as the satellite oscillates 
between pericenters and apocenters. But these short-time-scale variations are smoothed 
out and do not contribute to the evolution of the orbital angular momentum. In fact our 
results arc very insensitive to the number of pericentric passages for between a few to 
a few hundred. The disk and bulge are not important for our conclusion because we predict 
mainly the orbital decay in the outer halo where j > jdisk = 3000kpc km/s. Inside 15 kpc, 
our estimation of dynamical friction by an SIS model is inaccurate only for satellites on low 
inclination orbits. If satellites come in random inclinations, it is more common to find high 
inclination orbits, for which our models should be fairly accurate even inside 15 kpc. 

5.2. Orbital decay of the progenitors of the Sgr and the Canis Major dwairf 

galaxies 

The Sgr dwarf and the Canis Major dwarf are the closest known dwarf galaxies, about 
15 kpc from the center of the Milky Way and at the edge of the Milky Way disk. The Canis 
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Major is on a (direct or retrograde) orbit slightly inclined from the plane of the Milky Way, 
and the Sgr is on a nearly polar orbit. Both orbits have a fairly low angular momentum 
with 5* ~ 20 kpc; the data on Sgr are more complete, and show that it oscillates between 10 
kpc periccntcr and 50 kpc apocentcr. Both contain several globular clusters. It is possible 
that the two dwarfs are the stripped-down version of a more massive object, which has 
dynamically decayed from the outer halo. 

Interestingly a possible extension of the Sgr has been reported recently in the SDSS data 
near the position of the outer halo globular cluster NGC2419 (Newberg et al. 2003). There is 
a stream-like enhancement of halo A-colorcd stars at the SDSS magnitude of qq = 20.3 in the 
plane of the Sgr's orbit, corresponding to a distance of 90 kpc. If this is true, it would imply 
that the Sgr has changed its orbits in the past Hubble time. There are two possible ways 
that this could happen. One is that the Sgr's orbit has been deflected by a massive satellite, 
such as the LMC or SMC. Indeed the orbits of the Sgr and the Magellanic Clouds do overlap 
at the Galactic poles, and simple timing arguments show that these systems encounter or fly 
by each other about 2.5 Gyrs ago at about 50 kpc on the North Galactic Pole if the rotation 
curve of the Milky Way is nearly flat (Zhao 1998). The problem of this solution is that it is 
rare for the Sgr to receive a strong enough deflection to bring down its orbit. 

Another solution is that the Sgr has been a more massive system, which orbital decayed 
from the outer halo (Jiang & Binney 2000). This solution requires a large amount of massloss. 
Suppose that the progenitor of the Sgr has a mass comparable to the LMC, IO^^Mq, and 
was on an orbit with similar pericenter and apocentcr as the LMC, i.e., between 40 kpc and 
100 kpc, then it would nicely explain the stream near NGC2419 as debris of the Sgr at an 
earlier time, when the Sgr was on a larger orbit with apocenter greater than 90 kpc (cf. Fig. 

2). 

5.3. The orbit of a;-Centauri and its progenitor 

We have mainly concentrated on the problem of getting rid of a satellite's angular 
momentum if it starts with a high angular momentum or orbital size 5*0 S> 15kpc. What 
would be the remnant distribution if a satellite is born with an initial orbital size Si < 15kpc? 
The stars in such a system are assembled in the inner halo from the start, e.g., by colliding 
an inf ailing gas cloud with the protogalactic gas clouds in the inner halo, (Fellhauer & 
Kroupa 2002). Or the stars form from extragalactic gas and descend on a very radial orbit, 
penetrating the inner 15 kpc of the host halo from its very flrst pericentric passage. An 
intriguing example is a;-Centauri . Unfortunately our analytical model is not suited for this 
system because it is presently on a low-inclination eccentric retrograde orbits between 1 
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and 6 kpc from the Galactic center (Dinescu et al. 1999), so the contribution of dynamical 
friction by the disk is important. Also hydrodynamical friction with the disk gas can play 
a role for an early-on partially gaseous satellite. Nevertheless, if one applies simplisticly 
the tidal massloss and ChandraSckhar's dynamical friction in a spherical halo,i one finds 
that while it seems easy to peel off a satellite galaxy to make a central star cluster, most 
simulations produce remnants on much larger orbits than a;-Centauri (Zhao 2002). It seems 
some fine tuning is required to select progenitors on very low angular momentum and/or low 
energy orbits: the initial angular momentum needs to be low enough for the progenitor to 
penetrate into the inner halo or the present location of a;-Centauri on its very first pericentric 
passage. This means that the initial orbital size Sf of a;-Centauri is in between the present 
value of a;-Centauri Sq ~ 1.25kpc and the boundary of the inner halo it!disk = 15kpc, or 
mathematically 

1.25 kpc < < 15 kpc. (22) 

Most recently there have been several very encouraging attempts to model the dynam- 
ical and star formation history of a;-Centauri by nearly self-consistent N-body simulations 
(Mizutani et al. 2003, Tsuchiya et al. 2003, Bekki & Freeman 2003). All are able to produce 
a reasonable mass and orbit of the a;-Centauri simultaneously after some trial and error with 
the initial parameters of the progenitor; many initial conditions lead to remnants, unlike 
cJ-Centauri , beyond lOkpc of the Milky Way center. The favored initial orbit has a small 
orbital size Si. according to Tsuchiya et al. ji = 60kpc x 20kms~^ = 1200kpckms~^ (or 
Si — 6kpc) and according to Bekki & Freeman ji = 25kpc x 60kms~^ = 1500kpckms^^ (or 
Si — 7.5kpc). The small orbital size seems consistent with our expectation (cf. Fig.3). 

Tsuchiya et al. launch satellites with various initial mass (0.4 — 1.6) x 10^'^ M0 and with 
cither a King profile or a Hernquist profile from 60 kpc from the Milky Way center. They 
choose well-aimed nearly radial orbits, with an initial perigalactic radius about 1 kpc, much 
more radial than the present eccentric orbit. Massloss in their King model are similar to our 
exponential massloss models {n — 00): rapid in the beginning, and log(m) is roughly linear 
with time up to a mass of IO^Mq when the satellite has too httle mass to proceed with the 
orbital decay. Massloss in the Hernquist model is closer to a n = 0.3 model, linear in the 
beginning and rapid just before complete disruption. 

The progenitor in Tsuchiya et al.'s best simulation is a two-component "nucleuated" 
model with a sateUite of 0.8 x IO^^Mq with a Hernquist profile of half-mass radius 1.4kpc 
plus a nucleus modeled by an extended-particle of IO^Mq with a half mass radius of 35pc. 
A closer examination of this nucleus, however, reveals a subtle internally inconsistency in 
the initial condition. The nucleus is too fiuffy to be there in the first place. Inside the 
half-mass radius of 35 pc of the nucleus, the nucleus or extended particle has less enclosed 
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mass Mnuc('") than the Hernquist main body of the sateUite MHern('"), or mathematically, 

M„ud(r) = O.5xlO^M0 < Mnernir) = SxIO^Mq ( ——) = 2x10^ r = 35pc. 

\r + (v2 - 1) • 1414pc/ 

(23) 

So the "nucleus" is not a real self-gravitating identity, and should have been disintegrated by 
the tides of the Hernquist component if it was not kept as a rigid particle by construction. 
If this extended particle were removed from the simulation, there will be no cj-Centauri -like 
remnant left once the whole Hernquist satellite is disrupted by Galactic tides. One could 
stablise the initial condition by replacing the extended particle by a black hole or a dense 
cluster with half-radius less than 35pc, but then the size of a;-Centauri is not reproduced 
satisfactorily. In this aspect, the nucleated dwarf models of Bekki & Freeman are more 
self-consistent. 



6. Summciry 

We have explored the orbital decay for a dwarf satellite with a range of initial specific 
angular momentum and massloss history. We follow the evolution of satellites in the mass- 
distance plane, and find generally very little evolution of specific angular momentum by 
dynamical friction. The progenitors of inner halo globular clusters and substructures can not 
be born on orbits of comparable angular momentum as present-day halo satellite galaxies. In 
general, Gl-type of objects at relatively large radius (assume Gl is on a circular orbit) should 
be much more common than the cj-Centauri -type inner galaxy objects if tidal stripping of 
outer halo dwarf satellites has been important in the formation of massive galaxies. It is also 
difficult for a minor merger to bring a massive BH from outside the host galaxy all the way 
to the host galaxy center. Satellite remnants (BHs, globulars and streams) tend to hang up 
in the outer halo. 

I thank Ken Freeman, George Meylan, Jim Pringle, Floor van Leeuwen for enlightening 
discussions during the a;-Centauri conference at loA, and Oleg Gnedin, Mike Irwin, Pavel 
Kroupa, Mark Wilkinson for a careful reading of an earher draft. Special thanks to MianMian 
and YiYi (http://www.ast. cam. ac.uk~hsz/Twins/First Year. html) for not making excessive 
noise and comments during this work. 
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Fig. 1. — shows the satelhte mass as a function of time after rescahng by the initial mass 
TTii (from top to bottom n = 0.3, l.cxo), and the rotation speed Kir(^) {ioi top to bottom 
p = ^, |) for the past Hubble time after rescaling by the present rotation speed Vq. 
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Fig. 2. — shows the predicted evolution histories of a sateUite in the plane of its effective mass 
and characteristic orbital size S{t) (which is the specific angular momentum j(t) divided 
by a chacteristic velocity 200kms~^ irrespective of eccentricity of the orbit). Also indicated 
are the estimated orbital size and the mass of the satellite galaxies of the Milky Way and 
M31. The upper hatched horn-like area shows models with a massloss between exponential 
(upper dashed boundary, n = oo) and linear (lower solid boundary, n = 1); the potential 
grows moderately with p = \- The lower hatched horn- like area shows the predictions with 
a minor evolution of the potential (p = |) with a massloss between linear (upper dashed 
boundary, n — 1) and accelerated (lower solid boundary n — 0.3). A massive satellite starts 
at = 4 Gyrs from the upper right corner with an effective mass ^{ti) — rrii — IQ^^Mq with 
angular momentum ji = 250kms~^ x 50kpc, and ends with a mass of IO^Mq. For different 
assumptions of the massloss rate, the intermediate mass and position of the remnant are 
indicated with a time step of 1 Gyr. Note the failure to deliver remnants to the lower left 



corner. 
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Fig. 3. — shows the initial orbital size Si = ji/200 and the progenitor mass of a satellite 
to deliver a small remnant (lO^M©) to the inner galaxy (upper curves) and to the Galactic 
center (lower curves). An object in the upper left corner cannot evolve to a globular in the 
inner halo within a Hubble time. Different line types and shaded regions have the same 
meaning as in Figure 2. Also indicated are the present values for the most massive globular 
clusters of the Milky Way and M31 and a central million solar mass BH. 
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Fig. 4. — compares the circular velocity curves of observed dwarf galaxies (upper panel): 
DD0154 (Carnigan & Purton 1998), NGC3109 (Jobin & Carignan 1990), NGC5585 (Blais- 
Ouellette et al. 1999), NGC2976 (Simon et al. 2003) with those infered from our model of 
massloss and dynamical friction (lower panel). The dashed thick hne to the upper right is 
the Keplerian circular velocity curve of a point mass of IO^^Mq, to the lower left is sohd body 
curve with a crossing time 25 Myr. The three thinner curves are models with n — 0.3, 1, 30 
(up to down, marked in different symbols). All models start with an initial mass = 
IO^^Mq and on an eccentric orbit with peri-to-apo ratio 1 : 5 and ji = 50 kpc x 250kms~^, 
similar to the LMC. The host galaxy has a moderate growth of the potential well with p = ^■ 
All models have solid-body cores, each is labeled with the final orbital radius 5*0 of the fully 
disrupted satellite (30kpc, 45 kpc, 60kpc from up to down). 
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Fig. 5. — Similar to Fig.4b with the three shades zones being models with n = 0.3, 1, 30 and 
the upper and lower boundary of each zone being host halo models with p = ^ and P = |- 
All models end with an angular momentum jo — 15 kpc x 200kms~^ with pericenter outside 
the inner 15 kpc of the halo. 
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Table 1. Specific angular momentum (SAM) or orbital size of known satellites of the 

Local Group 



jo — rkpc X V km s ^ 


-^0 = 200 kms-^ (kPc) 


Object 


Ref 


5kpc X 50kms~^ 


1.25 


a;-Centauri , 


Dinescu et al. 1999 


16kpc X 260kms~^ 


20.8 


Sgr stream, 


Ibata et al. 1997 


60kpc X 250kms~^ 


75 


Magellanic stream, 


Kroupa & Bastian 1997 


70kpc X 200 km s"^ 


70 


Ursa Minor 


Schweizer et al. 1997 


138kpc X 310 km s"^ 


213 


Fornax 


Piatek et al. 2002 


150kpc X 20kms"^ 


15 


M31 stream'^ 


McConniche et al. 2003 


15kpc X 200 km s"^ 


15 


Canis Major dwarf galaxy 


Martin et al. 2003 



^Assuming NGC205 or M32 is near the pericenter of the stream orbit. 



